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F. ' 

.£h ! Proton-proton correlations and correlations of p-a, d-a, and t-a from spectator decays 

/\! following 197 Au+ 197 Au collisions at 1000 AMeV have been measured with an highly ef- 

ficient detector hodoscope. The constructed correlation functions indicate a moderate 
expansion and low breakup densities similar to assumptions made in statistical multifrag- 
mentation models. In agreement with a volume breakup rather short time scales were 
deduced employing directional cuts in proton-proton correlations. 
PACS numbers: 25.70.Pq, 21.65.+f, 25.70.Mn 

1. INTRODUCTION 

Understanding the simultaneous formation of many fragments helps to gain an insight 
into multifragmentation. Densities lower than ground state density of nuclei are a pre- 
requisite of statistical models describing multifragmentation |l]|| . While static statistical 
models form fragments in an expanded volume breakup, the dynamic statistical model 
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Figure 1. Proton-Proton correlation- functions for cen- 
tral collisions. The top panel shows directional selec- 
tions generated with harmonic weights, the bottom panel 
shows directional selection of pairs of protons, whose an- 
gular ranges of the angle between relative and sum mo- 
mentum are indicated in the figure. 



[[|,|J combines an expansion with surface emission until the volume breakup of the remain- 
ing source occurs. Interferometry-type methods are widely considered as a valuable tool 
in determining the space-time extent of such a source and recently, spectator remnants 
of the reaction Au + Au at 1 AGeV incident energy were found to break up at densities, 
considerable lower than ground state density |5j . In that analysis a instantaneous volume 
breakup was assumed. In such a breakup the mean time difference between the emission 
of two fragments is assumed to be short, due to fluctuations of the breakup time. 

In this article we report on results of correlation measurements of protons from the 
target spectator following the collisions of Au + Au at 1 AGeV incident energy. Coin- 
cidences between protons were used to construct correlation functions. The results are 
found to be consistent with low breakup densities with values close to those assumed in 
the statistical multifragmentation models and short emission time differences with values 
close to those anticipated for volume breakup. 



2. EXPERIMENT 

Targets, consisting of 25 mg/cm 2 of 197 Au were irradiated by an 1 AGeV Au beam 
delivered by the Schwerionen-Synchrotron (SIS) at GSI in Darmstadt. For the results, 
presented here, we employed one multi-detector hodoscope, consisting of a total of 96 Si- 
CsI(Tl) telescopes in closely packed geometry. The hodoscope covered an angular range 
Qiab from 122° to 156° with the aim of selectively detecting the products of the target- 
spectator decay. Each telescope consisted of a 300 /im Si detector with 30 x 30 mm 2 
active area, followed by a 6 cm long CsI(Tl) scintillator with photodiode readout. The 
distance to the target was 60 cm. The products of the projectile decay were measured 
with the time-of-flight wall of the ALADIN spectrometer || and the quantity Zbound 
was determined event-by-event. Zbound is denned as the sum of the atomic numbers 
Zi of all fragments with Zj> 2. It reflects the variation of the charge of the primary 
spectator system and serves as a measure for the impact parameter. 



3. DATA ANALYSIS 



3.1. Correlation functions 

The correlation functions were constructed dividing the spectrum of relative momenta 
of two coincident particles by the spectrum of pairs from different events. One observes 
at a relative momentum of q ~ 20 MeV/c a maximum of the correlation function, which 
is dominated by the attractive singlet scattering between two protons. The height of the 
maximum is inversely related to the relative distance between two protons and hence, the 
diameter of the source for simultaneous emission J7|||. Different emission times of the 
protons cause an elongation of the source in the direction of the hodoscope |J . Since two 
protons in the perpendicular direction are closer, they suffer due to the Pauli blocking a 
stronger suppression. We selected pairs of protons with directional cuts on the angle \& 
between the sum momentum P and the relative momentum q. The choice of the sum 
momentum as reference direction is motivated by the tilted source elongation with respect 
to the beam axis. The notion for longitudinal (P || q) and transversal (P _L q) correlation 
functions, therefore, differs from what is commonly used in high energy physics. Instead of 
using cuts on the angle \1/ between sum momentum and relative momentum we employed 
for the first time a harmonic weight for each event. We used cos 2 (\l/) and sin 2 (\l/) as 
weight for the longitudinal and transversal correlation function, respectively. 

The proton-proton correlation-functions are shown in Fig. [I] for the most central im- 
pact parameter bin. We used protons with energies E > 20 Mev in order to prevent 
distortions from sequential decay [ID]. The correlation functions show a pronounced peak 
due to the final state interaction at q = 20 Mev/c relative momentum. The longitudinal 
selection exhibits at lower relative momenta, q < 50 MeV/c, correlations which are ~ 5% 
stronger compared with the transversal selection. For a source with Radius R = 7 f m one 
can roughly estimate from the uncertainty relation the relative momenta where antisym- 
metrization effects should take place as q = v^l 2 " < v3 * ^/ r ~ 50 MeV/c. The bottom 
panel of Fig. [I] serves as comparison between the two methods for generating directional 
cuts. Here, pairs of protons were selected according to the orientation between sum 
momentum and relative momentum as it was done in analyses performed before ||. Em- 
ploying harmonic weights gives the advantage to use the full data set for one directional 
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Figure 2. Energy spectra of protons 
at angles of 130, 140, and 150 deg. in 
the laboratory system. 
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Figure 3. x 2 distributions as a function of radius and emission time of the emitting source for 
four impact parameter ranges indicated by the Z bound ranges in the panels. The lines in the 
left panel are explained in the text. 



selection, while for the hard cuts less than half of the events are used. In detector setups 
with limited angular acceptance a hard transversal cut limits the statistic for large rela- 
tive momenta and makes the normalization more difficult. This eg. can be seen from the 
larger error bars of the transversal correlation function in the bottom panel of Fig. |l| for 
relative momenta q > 60 MeV/c. Here, many of the proton pairs with large transversal 
relative momenta miss the hodoscope. 

3.2. Monte-Carlo simulations 

The analysis of the p-p correlation functions was performed with the Koonin-Pratt 
formalism J7|||. Particles were chosen to be emitted randomly from the volume of an 
uniform sphere and their velocities were sampled according to a Maxwell-Boltzmann dis- 
tribution. An additional velocity component was added, which simulates the Coulomb 
repulsion dependent on the location of the particle within the source. The strength of 
this component was chosen recursively, in order to agree with the later extracted source 
radii. From the measured proton-energy spectra (Fig. |2]) one can observe that the emis- 
sion probability of the source is non isotropic. We modeled this anisotropy by assuming 
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Table 1 

Source parameters for the MC-calculation of the protons: given are the temperature T, the 
Gaussian momentum width of the transversal movement of the source, the maximum momentum 
due to the Coulomb repulsion, the Gaussian width of the proton-multiplicity spectrum, and the 
velocity of the source in beam direction v s top- 



a sideward movement (bounce) of the source perpendicular to the beam, which causes a 
reaction plane. It influences slightly the correlation functions by the event mixing tech- 
nique which we used: the choice of a pair of protons from events with different reaction 
planes causes an additional small relative momentum between the momenta. Therefore, 
the relative momentum spectrum of the mixed events is slightly shifted to larger values 
compared to the spectrum of the coincident proton pairs. The correlation function then 
shows a slightly decreasing slope. Since our experiment lacks an azimuthal symmetric 
solid angle coverage, we were not able to determine the reaction plane and to rotate the 
different events in the same plane to avoid these effects [11]. We rather preferred to in- 



troduce the reaction plane in the Monte-Carlo calculations. In order to simulate the hit 
probability of the hodoscope for different bounce values, we adjusted the multiplicity of 
the event such that it resembles the measured proton-multiplicity spectrum after applying 
the experimental filter. The extracted source parameters are given in Tab. [I]. 



4. RESULTS 

We varied the radius of an uniform density distribution and the Gaussian emission 
time of the protons. The time distribution is guided by the idea of a statistically in- 
dependent emission of different emitters within the source, rather than the decay of an 
excited state (see eg. [0). The simulated correlation functions were used to perform a 
X 2 -test algorithm in the relative momentum region of 10 < q < 35 MeV/c. The results 
are presented in Fig. |||. The grey scale of the shading denotes within which standard 
deviation the space-time point is. The minima in Fig. ^ show approximately constant 
source radii and surprisingly short emission times. For the most peripheral bin a radius 
and emission time could due to the low statistics of the correlation function not be de- 
duced. The experimental correlation functions (symbols) and the simulations (lines) for 
the best x 2 are shown in Fig. f|, left panel. They agree well with their simulations. The 
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Figure 4. The left panel 
shows experimental longi- 
tudinal (open symbols) and 
transversal (closed symbols) 
correlation functions and 
results of MC-simulations 
(lines). The right panel com- 
pares the their experimental 
ratios (symbols) with the 
results of the MC-simulation 
(lines). 
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minimum x 2 values per degree of freedom are within the range of 1.3 < x 2 /d.o.f. < 2.4. 
The right panel in Fig. [| presents the ratios between longitudinal and transversal corre- 
lation functions. One recognizes the weak enhancement of about 5% above unity (dashed 
line) of data (symbols) and simulations (solid line) for relative momenta q < 50 MeV/c. 
For low values, q < 20 MeV/c, this enhancement becomes suppressed by the Coulomb 
repulsion between the protons. Although the sensitivity of the directional cuts extends 
up to relative momenta of q < 50 MeV/c, we restricted our fit region to lower values 
of 10 < q < 35 MeV/c. In this region the shape of the simulated correlation functions 
describes well the data. The ratios in the right panel of Fig. ^ agree, however, similar 
well for the larger relative momenta. Using a fit region of 10 < q < 50 MeV/c yields radii 
within AR = 0.2 fm and emission times within At = 2 fm/c to those results extracted 
from the narrower relative momentum region. 

5. DISCUSSION 



The last collision of the two particles within the surrounding nuclear debris determines 
their freeze-out radius and is believed to occur when their mean free path gets larger 
than the system itself: the system becomes below a certain density transparent. We used 
reconstructed residue masses |T3| to calculate the freeze out densities (Fig. g). They are 
lower than the ground-state density of nuclei and therefore, show expansion. The strong 
variation of the freeze-out density is caused by the variation of the mass of the spectator. 
The freeze-out densities shown as open symbols were analysed in a previous work || with 
the assumption that the particles are emitted instantaneously. Correcting the source size 
with help of the directional cuts for the emission times of the protons results in densities 
(closed symbols) which are only slightly larger compared to the zero time assumption. The 
expected higher freeze-out densities were lowered by two effects, the collective movement 
of the source and by the initial Coulomb correlations within the source. In the left 
panel of Fig.[| we show for 20 < Zbound < 40 the simulation of a purely thermal source 
(thin lines) with the same radius. Guided by the freeze-out conditions of a statistical 
multifragmentation model [ffl] we tried to understand the variation of the freeze-out density 
with the impact parameter: assuming a constant crack width between the fragments, the 
freeze-out density becomes dependent on the multiplicity M. Hereby, the diameter of a 
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Figure 5. Freeze-out densities from p- 

p, p-a, and d-a correlation functions (open 
symbols) were analysed assuming an instan- 
taneous emission of the particles ||. The 
closed symbols represent the freeze-out densi- 
ties extracted from p-p correlation-functions 
assuming a Gaussian time distribution. The 
lines are explained in the text. 



nucleus in one dimension is compared to the diameter of an expanded nucleus with M 1//3 — 1 
cracks with the width 2d. This one dimensional freeze-out radius is independent on the 
fragment size and only an approximation in three dimensions. In this simplified picture 
the freeze-out volume over the ground state volume is Vf reeze /V = (R +d(M 1 ^ 3 — 1)) 3 /Rq. 
We used the event multiplicity from Ref. ]n| for different species of fragments. For the 
calculation of the primordial multiplicities we assume a simple sequential decay estimate: 
we observe in energy spectra for neutrons, protons, and alpha-particles a low energy 



component, which for the protons scales with the biggest charge in an event [15| and is 



likely due to sequential decays. We, therefore, select for those particles the high energy 



component of the spectra. For the crack width we assume two extreme assumptions Jy 
for a gas of nucleons the crack width 2d is in the order of 1.4 fm, the range of nuclear 
interaction. For an event with heavier fragments the fragments surfaces are separated by 
the distance where the nuclear attraction equals the Coulomb repulsion (2d = 2.8 fm). 
Both assumptions are plotted in Fig. |5] as dashed lines. They exhibit the expected 
decrease of the freeze-out density with centrality. The extracted proton-proton freeze-out 
densities corrected for emission times follow closely within 20% the predictions for a crack 
width of 2d = 1.4 fm. 

The other scattering states beside of the most peripheral reactions are in between 
of both predictions. These scattering states are, however, not corrected for emission 
times. The radius R T =o with inclusion of the emission times r can be estimated from 



the uncorrected radius R e jf = JR^ =0 + (vt) 2 , with the thermal velocity v for a fragment 



with temperature of T = 17 MeV [|15[ . This estimate is plotted as solid and dashed lines 
in Fig. |3] for Zbound < 20 for radii from pa- and da-correlation functions, respectively. 
Emission times between r = 10 — 30 fm/c are needed to get a freeze-out radius comparable 
with that extracted from the emission-time corrected pp-analysis. 

The times between the emission of protons are in the order of the passing time of the 
projectile through the target. We, therefore, cannot exclude that the protons come from 
the fireball, which grows with the centrality of the reaction. Attempts to fit the energy 
spectra in Fig. |] with an fireball localized in momentum space failed. Nonetheless, the 
intriguing similarity of the impact-parameter dependence of the heavier scattering states 
(pa, da) which differ only by a common factor from the pp freeze-out densities indicate 
the emission from the target spectator. 

It is interesting to point out that the average expansion from the RMS radius at 
ground-state to the freeze-out density within our emission times is much too high. It 
is (/?) = 0.19 — 0.29 c for the most central bin and disagrees with the observed energy 
spectra. A low expansion velocity of about E/A = 0.5 MeV (0.03 c), suggested by other 



studies [[nj, has the consequence that the nucleus first expands for t w 100 fm/c and then 
decays within the extracted short times of r = 10 — 20 fm/c. A reaction where the spec- 
tators are diluting each other by scattering out the fireball nucleons would explain the 
large radii without substantial amount of expansion. Calculating the expansion velocity 
of a nucleus starting to expand from the ground-state radius of a 197 Au nucleus yields an 
average expansion velocity of (/?) = 0.04 — 0.09 c for the most central bin. 



6. CONCLUSION 

We constructed correlation functions for pairs of protons detected at backward angles 
in the reaction Au+Au at 1000 AMeV incident beam energy. Using high energy cuts of 



E > 20 MeV we selected protons which are only little affected by sequential feeding [|1C 
Comparing the results with Monte-Carlo simulations within the Koonin-Pratt formalism 
fairly constant freeze-out radii of R ~ 8 fm are deduced and emission times of r = 10 — 15 
fm/c are surprisingly short. The extracted radii are larger than the ground state radii 
of target spectators and show expansion. The mass dependence of the decaying target 
spectator with impact parameter lead to non constant freeze-out densities which decrease 



with centrality. Employing a multiplicity dependent freeze-out picture |TJj we are able 
to describe the absolute value of the freeze-out densities as well as the impact parameter 
dependence. Because of the short emission times in the order of the passing time of both 
spectators we cannot exclude that the protons come from first stage scattering of the 
nuclear cascade. However the impact parameter dependence of the freeze-out densities 
of heavier scattering states from a former analysis PJlOl indicates the emission from the 
target spectator. A more detailed analysis of these heavier scattering states is necessary 
to understand their formation and freeze-out properties. 
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